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There are several approximations to the exchange-correlation functional in density-functional the¬ 
ory that accurately predict total energy-related properties of many-electron systems, such as bind¬ 
ing energies, bond lengths, and crystal structures. Other approximations are designed to describe 
potential-related processes, such as charge transfer and photoemission. However, the development of 
a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to 
address it, we employ in the current work the ensemble generalization procedure proposed in Phys. 

Rev. Lett. 110 , 126403 (2013). Focusing on the prediction of the ionization potential via the highest 
occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the 
local spin-density approximation, semi-local generalized gradient approximations, and global and lo¬ 
cal hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. 

We hnd that the aforementioned ensemble generalization systematically improves the prediction of 
the ionization potential, for various systems and exchange-correlation functionals, without compro¬ 
mising the accuracy of total energy-related properties. We specifically examine hybrid functionals. 

These depend on a parameter controlling the ratio of semi-local to non-local functional components. 

The ionization potential obtained with ensemble-generalized functionals is found to depend only 
weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus 
eliminating one aspect of the so-called ‘p^-rarneter dilemma’ of hybrid functionals. 


I. INTRODUCTION 


Modern density-functional theory (DFT), based on the 
theoretical foundation laid by Hohenberg, Kohn, and 
Sham 0,0 in the 1960s, in principle provides an ex¬ 
act framework for treating the many-electron problem. 
Within this framework, the electron-electron interaction 
is expressed via the exchange-correlation (xc) energy 
term, Excjn], which is a functional of the electron den¬ 
sity n(r) 0-0. In practice, one has to approximate this 
energy contribution, aiming at a numerically efficient yet 
accurate description of the electronic structure of var¬ 
ious many-electron systems, such as atoms, molecules, 
and solids. 

During the eventful history of DFT, many density 
functional approximations (DFAs) to the exact xc energy 
were developed 00 . The performance of each of them 
can be evaluated from a two-fold perspective: On the one 
hand, ground-state quantities such as binding energies, 
bond lengths, and crystal structures are closely related 
to the total energy of the system and hence to the ap¬ 
proximate xc energy Exc itself. On the other hand, there 
also exists a range of physical properties and processes, 
whose description is substantially influenced by the xc po¬ 
tential, nxc(r) := ^£’xc[n.]/5n(r). Prominent examples for 
the latter category are charge-transfer and ionization pro¬ 
cesses, as well as the description of photoemission spec¬ 
tra. Here, especially the approximate interpretation of 
Kohn-Sham (KS) eigenvalues as a physically meaningful 
density of states lies at the focal point of ongoing research 
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(see Ref. [§| and references therein). 

For applications to real materials, one would wish to 
have a DFA with a good performance from both perspec¬ 
tives. However, this is not the case for many existing 
DFAs. The development of such a DFA is a long-standing 
challenge, as discussed below. 


One of the exact relations in KS-DFT is the ioniza¬ 
tion potential (IP) theorem, —£ho = I UM , which re¬ 
lates the highest occupied (ho) KS eigenvalue, £ho, to 
the IP, I, i.e., the removal energy of one electron from 
an A"-electron system. Consequently, there are two fun¬ 
damentally different ways to obtain the IP in DFT: by 
evaluating —Eho or (with more computational effort) by 
calculating the energy difference of the ionized and neu¬ 
tral system, with A^q — 1 and Nq electrons, respectively: 
I^ASCF = E{No — 1) — E{No). This is usually referred to 
as the ASCF approach. Note that while the latter relies 
on accurate total energy values, the former relies on the 
potential to yield an accurate ho energy level. 

The ASCF approach, which has been extensively used 
since the early days of DFT [311, usually yields an 
IP with a satisfactory accuracy of a few percent with 
respect to experiment, for atoms and small molecules, 
even with standard (semi-)local DFAs, such as the lo¬ 
cal spin-density approximation (LSDA) [2g or the gen¬ 
eralized gradient approximation (GCA) [27|. Compari¬ 
son of —£ho to the experimental IP, however, shows poor 
correspondence for many DFAs. For example, with the 
aforementioned (semi-)local approximations, one can ob¬ 
serve an underestimation of up to 50%, as manifested in 


e.g. 


This failure has been related to various systematic 
shortcomings of existing functionals. First, the one- 
electron self-interaction problem, i.e., the fact that for 
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many DFAs the Hartree energy Ah is not canceled by Avc 
if evaluated on one-electron ground-state densities j22j |. 
is well known to have a large impact on the quality of 
KS eigenvalues 0. Second, many DFAs show a potential 
that features an incorrect long-range asymptotic behav¬ 
ior, again with negative effect on the interpretation of 
KS eigenvalues (see, e.g., and references therein). 
Note that while the two issues are related in a physical 
sense, they are not the same and their connection in the 
construction of reasonable DFAs is far less obvious [s^. 
Another shortcoming affecting the KS eigenvalues is the 
deviation of the total energy curve, E{N), as a function of 
the number of electrons, N, from piecewise linearity, for 
fractional N (see, e.g., [14 [35l - l4§ | 1. In the literature, this 
phenomenon is sometimes referred to as many-electron 
self-interaction [ssl - fs^ or as a (de-)localization er¬ 
ror [H, [dlj, |4^ . For (semi-)local functionals, one obtains 
a convex energy curve rather than a straight line. Such a 
deviation reflects negatively also on systems with an inte¬ 
ger N: it leads to disagreement between the IPs predicted 
by the ASCF method and those predicted via —Sho- This 
happens because the slope of the energy curve (to the 
left) equals, according to Janak’s theorem [d^, the ho KS 
eigenvalue. Therefore, even when one is interested only 
in closed systems, with an integer number of electrons, it 
is important to tackle the problem of lack of piecewise- 
linearity in order to obtain a physically meaningful value 
for Eho- 

There exist many approaches to address the aforemen¬ 
tioned shortcomings and obtain accurate results for the 
IP via Eho- Self-interaction correction [2^ schemes lead 
to a significant improvement in the interpretation of KS 
eigenvalues [H^. Yet, their performance for ground-state 
energetics is debatable [Sll - I^ . Approaches that approx¬ 
imate directly the xc potential [56l - l59l | yield eigenvalues 
that satisfactorily reproduce the experimental IP, due to 
modihed long-range properties of the potential. However, 
for these functionals total-energy related quantities are 
not accessible [hOl - l^ . New types of GGAs can yield sig¬ 
nificantly improved potential properties [63l465l |. but at 
the cost of being less accurate for total energies. Global 
hybrid functionals (66l - [7ll| , which linearly combine (semi- 
)local xc energy components, with a weight (1 — a), and 
exact exchange (EXX, i.e., the Fock-integral evaluated 
with KS orbitals), with a weight a, mitigate the one- 
electron self-interaction error and often yield an excel¬ 
lent description of properties related to the total energy, 
for a Ri 0.25. However, since the self-interaction is only 
partly canceled, the KS potential falls off too quickly in 
the asymptotic limit, and —£ho is typically far from de¬ 
scribing experimental IPs. 

Nevertheless, with global hybrids it is possible to find 
a value of a such that the global hybrid will produce a 
piecewise linear energy curve, and therefore an improved 
value for —Sho- This happens because for fractional N the 
non-local EXX component of the glob al hybrid produces 
a concave energy curve (see, e.g. [37|), while the (semi- 
)local components usually cause a convex energy curve. 


which therefore cancel each other. However, this cancel¬ 
lation is achieved with values of a « 0.75 [72 - l74l |. which 
in most cases significantly compromises t he p erformance 
of the functional for other quantities [H, . 

This creates what we call the ‘parameter dilemma’: 
while an accurate description of energy-related quanti¬ 
ties requires a certain value for the functional’s parame¬ 
ter, the accurate description of potential-related quanti¬ 
ties requires a different value, and there is no value that 
provides a satisfactory description of both [78l - [80| . 

Local hybrid functionals (sil - l^ aim at preserving the 
good energetics of global hybrid functionals while reduc¬ 
ing the self-interaction error by introduction of a more 
flexible, space-dependent mixing of (semi-)local and non¬ 
local components (see [s^ for an overview and discus¬ 
sion). However, we recently illustrated using a specially 
constructed local hybrid functional, termed ISOcc in the 
following, that the aforementioned ‘parameter dilemma’ 
persists also there [s^ . Similarly, in range-separated hy¬ 
brids (RSHs) the values of the range-separation param¬ 
eter have to be different to accurately reproduce e.g. at¬ 
omization energies and ionization potentials [^ . 

It has been recently shown [s^ [sj that an alterna¬ 
tive way to improve the prediction of the IP via —Ehn is 
given by employment of the ensemble approach [3, ISSMOOl l 
in KS-DFT. This approach allows for the generalization 
of the Hartree and xc functionals for fractional N such 
that the piecewise linearity behavior of the total energy is 
restored, to a large extent. As a result, better correspon¬ 
dence of £ho to the experimental IP and to the ASGF 
value is achieved, as demonstrated for the H 2 molecule 
and the G atom with the LSDA. 

Here we emplo y th e ensemble generalization procedure 
proposed in Ref. |87| to the Hartree and common approx¬ 
imate xc functionals, aiming to address the aforemen¬ 
tioned challenge of simultaneous prediction of energy- 
related and potential-related properties with one DFA. 
Focusing on the prediction of the IP via the ho KS eigen¬ 
value, we examine a variety of xc approximations: the 
local spin-density approximation, semi-local generalized 
gradient approximations as well as global and local hy¬ 
brids. Results for a representative test set of 26 light 
diatomic molecules and single atoms are presented. We 
find that the ensemble generalization systematically im¬ 
proves the prediction of the IP, for a wide variety of sys¬ 
tems and xc functionals, changing the general tendency 
from under- to a small overestimation, compared to ex¬ 
periment. This improvement is achieved without any 
change in total energy-related properties. For hybrids 
that include a parameter, the IP obtained with ensemble- 
generalized functionals is found to be only weakly de¬ 
pendent on the parameter value, contrary to common 
experience with non-generalized hybrids. Thus, the en¬ 
semble approach eliminates one aspect of the ‘parameter 
dilemma’. 
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II. THEORETICAL BACKGROUND 

For completeness, we briefly present the ensemble 
generalization to the approximate Hartree-exchange- 
correlation (Hxc) density functional, focusing on its in¬ 
fluence on the highest occupied KS energy level, £iin- A 
complete derivation can be found in Refs. 1^. 

First, we formally consider a system with a fractional 
number of electrons, N = Nq — 1 + a, where Aq S N and 
a € [0,1], so that a = 1 corresponds to a neutral and 
a = 0 to a singly-ionized system. Subsequently, we take 
the limit 0 ^- 1 “, focusing on neutral systems with an 
integer number of electrons. 

In a landmark article, Perdew et al. have shown 
that the zero-temperature ground state of an interacting 
many-electron system possessing a fractional N should 
be described by an ensemble state @. This state is a 
linear combination of the pure ground states for Nq — 1 
and Nq electrons, with the class i cal s tatistical weights of 
(1 — a) and a, resDectivelv [l34l | . |l35l |. The ground-state 
energy of this ensemble state has then been shown to be 
equal to E{N) = (1 — a)E{No — 1) -f aE{No), i.e., it is a 
piecewise-linear function of iV Q . This result is a general 
one, applying to any many-electron system. Therefore, 
in principle it trivially carries over to DFT, because if the 
exact exchange-correlation functional is used, DFT-based 
energies must reproduce the all-electron ones. 

As mentioned in the introduction, in practice approx¬ 
imate density functionals often exhibit significant devi¬ 
ations when trying to describe a quantum system with 
fractional N in KS-DFT, and it has been traditionally 
assumed that this is just another manifestation of the 
approximate nature of the functional used. However, in 
Ref. [13 it was pointed out that much of this deviation is 
due to the fact that the pure-state exchange-correlation 
expression is used for both integer and fractional densi¬ 
ties, whereas the fractional KS system must itself be in an 
ensemble state. This happens because the number of par¬ 
ticles in the KS system equals the number of electrons in 
the real, interacting system, and is also fractional. There¬ 
fore, the ground state of the KS system must also be ex¬ 
pressed as an ensemble of the {Nq — I)- and Nq- KS states, 
obtained from the same KS potential, even if one uses 
an approximate functional. Ref. H^j therefore suggested 
that any approximate Hartree-exchange-correlation func¬ 
tional can be generalized for an ensemble ground state us¬ 
ing ensemble state theory ffo r other recent uses 

of the ensemble approach see Refs. |46l . IMlIil)- Perform¬ 
ing an ensemble average of the many-electron Coulomb 
operator W = ^ I''* system, it 

has been found that the pure-state Hxc energy functional 
can be generalized to ensemble states in the following 
form: 

Ae.Hxc[n(“)] = (1 - a)KHxc[pL“^] + aEn.c[pi% (1) 

which is exact for the Hartree and exchange components 
and approximate for the correlation. Here, the index 
e- indicates that the functional is ensemble-generalized, 


Ahxc is the pure-state Hxc functional, pp“^(r) is defined 
as the sum of the first Nq p KS orbitals squared: 
Pp“^(r) = where p = -1 or 0, and 

n(“)(r) = (1 — (r) -|- apQ°'\r) is the ensemble-state 

electron density. When N is an integer, i.e., a assumes 
the value of 0 or 1, the Hxc energy reduces to that ob¬ 
tained from the underlying pure-state Hxc functional. 
Therefore, ensemble-generalization does not affect the to¬ 
tal energy at integer N. The generalization in Eq. 0 is 
applicable to any xc functional and makes the Hartree 
and the xc energy components explicitly linear in a. How¬ 
ever, there may still remain an implicit non-linear depen¬ 
dence of -E'e_Hxc[’^^“^] on a, because the KS orbitals them¬ 
selves, (p-“^(r), and consequently Pp°‘\r) and i?Hxc[Pp°'^], 
may depend on a. The dependence of on a arises 

from the fact that the KS orbitals are expected to relax 
as one varies a from 0 (positive ion) to 1 (neutral system) 

[illlllil. 

Importantly, Eq. © is derived by considering the 
generalization of pure-state functionals to ensemble 
states, without assuming anything a priori about 
piecewise-linearity, because it also applies to approxi¬ 
mate exchange-correlation functionals. Nevertheless, in 
Refs. [13, 113 it has been shown that by employing 
Eq. ([T]) the energy curve E{N) satisfies the piecewise- 
linearity criterion much more closely, being slightly con¬ 
cave. The concavity is related to the above mentioned 
implicit non-linear dependence of the energy on a. An¬ 
other perspective on this approximate piecewise-linearity 
can be obtained from the fact that Eq. o can be de¬ 
rived, with some further approximations, from different 
schemes that attempt to enforce piecewise-linearity ex¬ 
plicitly B B113 Ellll Hi ■ 

Due to the fact that the slope of E{N) changes for all 
a, iy lading a —>• I“, it follows from Janak’s theorem [Hj ] 
[l36j| . which identifies dE/dN with £ho, that the ho en¬ 
ergy level has to change, too, even for a system with an 
integer N . This change is obtained in practice from an 
ensemble generalization of the KS potential, as explained 
below. 

The KS potential is expressed as Ue-Ks(r) = Uext(r) + 
?;e-HxcM(r), where Uext(r) is the external potential and 
?^e-HxcM(r) := ^Ae-Hxc/<^’^(r) is the ensemble-generalized 
Hxc potential. Eor the limit a —>■ 1“, this potential re¬ 
duces to a sum of two terms: Ue-HxcH(r) = uhxc M (r) -f 
uo[n] “ the usual pure-state Hxc potential, UHxc[?^](r), 
and a spatially uniform term, uo[n], which can be writ¬ 
ten as 

VQ[n] = EhxcM - KHxc[n - 

- j |v5ho(I■)P^^Hxc[n](I•)d^^. (2) 

Here and below the superscript (a) is dropped at the limit 
a ^ 1“ for brevity. Note that the ensemble-generalized 
KS potential does not vanish at r —>■ oo, but asymptot¬ 
ically approaches Uo[n]. We stress that Uo[n] is a well- 
defined, rather than arbitrary, potential shift. It must 
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be taken into account for the ensemble-generalized func¬ 
tional in order for the ho KS eigenvalue to equal dE/dN, 
i.e., to obey Janak’s theorem. Note that the shift dis¬ 
cussed here is different from the one recently proposed 
by Zahariev and Levy (9^. As clarified in Ref. [o^, in 
Ref. the potential shift makes the energy of the KS 
system equal the energy of the interacting system. Here, 
however, the shift emerges naturally from the ensemble 
treatment and is essential to obtaining results that are 
consistent with Janak’s theorem. Also note that while 
the result above has been presented in a spin-independent 
form for simplicity, in practice, in spin-dependent calcu¬ 
lations, there exist potential shifts Vq to both spin chan¬ 
nels cr =t, 4^- Calculating Vq with Eq. we take the 
ho level to be the highest occupied level in the cr-channel 
considered (noted as cr—ho). In the following, however, if 
not stated explicitly otherwise, when mentioning the ho 
level we refer to the global ho: Eho = max^ e^o, i.e. the 
one of the two cr-ho levels which is higher in energy; the 
same applies for the ensemble-generalized ho level, £e-ho- 
To summarize, as a result of the approximate ensem¬ 
ble generalization of the Hxc functional (Eq. ([T])) 
in the limit of integer N the KS potentials exhibit spa¬ 
tially uniform shifts Vq, such that all KS eigenvalues of 
the same spin channel are shifted by the same value (see, 
e.g.. Fig. 4 in Ref. 93|. The KS orbitals, and as a result 
the density and the total energy, are not changed and 
remain the same as those obtained with the underlying 
Hxc functional. Furthermore, because all eigenvalues are 
shifted by the same amount, eigenvalue differences (as 
well as quantities based on them, e.g., in linear response 
time-dependent DFT [9^) are not affected either. There¬ 
fore, the cr—ho energy levels of the ensemble-generalized 
functional can be expressed as + Vq, being 

a sum of the cr-ho level that emerges from a standard 
KS-DFT calculation prior to the ensemble generalization 
and the potential shift of the relevant spin channel, cal¬ 
culated according to Eq. ([2]). Comparing both Ee-ho and 
£ho to experimental IPs and —/ascf is the main subject 
of Sec. [IS 


III. COMPUTATIONAL DETAILS 


We concentrate on a relatively elementary, yet chemi¬ 
cally representative, set of systems, consisting of 18 light 
diatomic molecules: H 2 , LiH, Li 2 , LiF, BeH, BH, BO, 
BE, CH, CN, CO, NH, N 2 , NO, OH, O 2 , EH, F 2 , and 
their 8 constituent atoms. The simplicity of the systems 
allows us to keep computational costs low and to refrain 
from introducing additional sources of error, e.g., search¬ 
ing for an optimal geometry in systems with many de¬ 
grees of freedom. At the same time, systems of single-, 
double-, and triple-bond molecules as well as atoms (no 
bonding) are included in the test set, which makes the 
set representative of more complicated systems, as shown 
in previous work (see, e.g., [^, js^)- 

All calculations were performed using the program 


package DARSEC [13, an all-electron code that al¬ 
lows for electronic structure calculations of single atoms 
or diatomic molecules on a real-space grid represented 
by prolate-spheroidal coordinates. DARSEC allows one to 
solve the KS equations self-consistently for density- as 
well as orbital-dependent functionals. For the latter, a 
local, multiplicative xc potential is obtained by employ¬ 
ing the KLI a pproxima tion to the optimized effective 
potential (OEP) [l00l - ll02 | formalism. Us e of this ap¬ 
proximation has been justified in Ref. [I03j| for the EXX 
functional and in Ref. [s^ for the ISOcc local hybrid 
functional. 

For all systems, an accuracy of 0.0005 Ha in the total 
energy and in the ho KS eigenvalue has been achieved 
by appropriately choosing the parameters of the real- 
space grid and by iterating the self-consistent DFT cycle. 
For m olecules t he bond length was taken from experi¬ 
ment [lol flo^ . Differences due t o at omic relaxation 
were found to be insignificant |l37j| . |l38j . The net spin 
of the neutral systems was also taken to be as in experi¬ 
ment. The spin configurations of cations (used below for 
calculating ionization potentials from total energy differ¬ 
ences) was obtained by removing an electron from the 
highest occupied orbital of the neutral. 


IV. RESULTS 

A. Effect of the ensemble correction - O2 as a 
prototypical case 

Previous work [13, US] has already demonstrated that 
the ensemble generalization of Eq. o significantly re¬ 
duces the deviation from the piecewise-linearity condi¬ 
tion for the total energy, i.e., greatly diminishes the de- 
localization error, and as a consequence eliminates the 
fractional dissociation error in diatomic molecules. Here 
we focus on the potential shifts (Eq. ([2|)) that emerge 
from the ensemble-generalization and their effect on the 
Kohn-Sham energy levels. In particular, we consider the 
prediction of the IP via £e-ho- 

For a clear understanding of the results presented in 
this paper, it is of advantage to first illustrate the ef¬ 
fect of the potential shift mechanism, given by Eq. 
on the eigenvalue structure of a particular system with 
an integer number of electrons. Here we provide a 
detailed presentation of a selected system - the O 2 
molecule, computed with the Perdew-Burke-Ernzerhof 
(PBE) GGA at its experimental bond length of 
2.2819 bohr. 

Due to its electronic ground-state configuration, 
this system must be treated in a spin-polarized formal¬ 
ism. Gonsequently, it provides an interesting example for 
how eigenvalues belonging to different spin channels are 
shifted when the corresponding ensemble potential shift, 
Vq, is applied. 

For this purpose, the positions of the highest occupied 
(^ho) lowest unoccupied (£[(,) KS eigenvalues for both 
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spin channels are depicted in Fig. [TJ The eigenvalues 
changed by the respective potential shift, i.e., = 

ego + Vq and eg_io = ejg + tig, as well as the negative 
of the experimental IP, —/exp, are also included in the 
figure. 


0.00 
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^e-ho_ _ 

^ ^ e-ho 


up channel down channel 


FIG. 1: Diagram of the highest occupied and lowest unoccu¬ 
pied KS-PBE eigenvalues of the O 2 molecule, for both spin 
channels, before and after applying the potential shifts of 
Eq. @, along with the negative of the experimental IP. All 
values are in Ha. 


It can be readily observed that the unshifted highest 
occupied eigenvalue of the up channel, = —0.251 Ha, 
which lies higher than its spin down counterpart, poorly 
reproduces the negative of the experimental IP of the O 2 
molecule. In fact, with PBE it und erestimates the exper¬ 
imental IP of /exp = 0.453 Ha [l05j| by 45 %, a value that 
is quite typical for other systems as well. However, after 
application of the potential shift, the highest occupied 
eigenvalue is = —0.526 Ha, i.e., the experimental 
IP is now overestimated by 16 %. As shown below, this 
is a typical result also for other systems and other func¬ 
tionals. 

From the results presented for the ensemble generalized 
PBE functional, the question of how other DFAs perform 
for the same system naturally arises. In particular, the 
change of the eigenvalues obtained with functionals con¬ 
taining a varying amount of non-local EXX is of great 
interest, as we know that a greater percentage of EXX 
already leads to a more accurate description of IPs via 
the highest occupied eigenvalue. 

Fig. [2] provides a comparison of the unshifted and 
shifted ho eigenvalues to the experimental IP for 02- 
Besides (se mi-~)local f unctionals such as the LSDA, PBE 
and BLYP we also ensemble-generalized the 

global hybrid functionals B3LYP [6^ and PBEh(a) [7flj | 
(employed within the KS scheme using the KLI approxi¬ 
mation), with a denoting the fixed amount of EXX com¬ 
bined with (1 — a) of PBE exchange and with full PBE 
correlation. 

The unshifted eigenvalues for the three purely (semi- 
)local functionals (LSDA, PBE, BLYP) underestimate 
the IP by ~ 45 %. After ensemble generalization, we 
observe an overestimation by ^ 16 %. It is instructive 
to check to which extent this overestimation comes from 



FIG. 2: Comparison of —£ho and —Ee-ho to the experimental 
IP of the O 2 molecule, calculated with different DFAs. The 
corresponding labels provide the relative deviation in percent. 


errors in calculating total energy differences that are in¬ 
herent to the underlying functional, and to which e xten t 
they arise from the ensemble generalization process [l39j| . 
Therefore, we compare the shifted and unshifted eigen¬ 
values to — /ascf- For O 2 computed with PBE, one ob¬ 
tains /ascf = 0.464 Ha, which deviates from experiment 
by only 2.3 %. With respect to this quantity, the un¬ 
shifted eigenvalue yields an underestimate of 46 %, while 
the shifted value overestimates it by 14 %. We there¬ 
fore realize that most of the discrepancy comes from the 
concavity that remains in the E(N) curve even after en¬ 
semble generalization. 


In the global hybrid, PBEh(a), increasing the intrin¬ 
sic amount of EXX signihcantly improves the correspon¬ 
dence of the unshifted eigenvalue to experiment. Due 
to the fact that the KS potential decays asymptotically 
more slowly with a growing value of a, the IP via —Eho 
is very sensitive to non-local functional components in¬ 
cluded. Changing from under- to overestimation, an op¬ 
timal description of /exp is reached for this system with 
a ~ 0.6. However, for ensemble generalized DFAs the 
value of —Ee-ho systematically overestimates the IP with 
respect to experiment for the O 2 molecule, regardless of 
the value of a, while at the same time being far less sen¬ 
sitive to the amount of non-locality in the functional ex¬ 
pression. While for ’’plain” PBE the relative error now 
reads ~ -1-16 %, it increases to -1-29 % when full non-local 
exchange combined with PBE correlation is used. The 
reason for this reduced sensitivity lies in the following 
mechanism: while the absolute value of £ho grows with 
increasing a, the potential shift vq is reduced, roughly 
commensurately, because the Hartree-|-EXX functional 
has zero potential shift [s^- 
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B. Evaluating the test set - a systematic study 


Following the illustration of the mechanism of the po¬ 
tential shift for a single system, we now focus on the 
mean discrepancy in the evaluation of the experimental 
IP via shifted and unshifted KS eigenvalues, for a variety 
of functionals (see supplemental material at [URL will 
be inserted by AIP] for detailed numerical data in tabu¬ 
lar form). We use the test set of systems introduced in 
Sec.|nl]as a basis for averaging. 

We emphasize that the eigenvalues shift is expected 
to improve the correspondence between the negative of 
the ho eigenvalue and the ionization energy obtained via 
the ASCF method, for a given DFA. We compare the 
shifted and unshifted eigenvalues to experiment, and not 
to ASCF values, relying on the aforementioned fact that 
the ASCF reliably describes systems of our test set, with 
small average relative errors: 3.4 % for PBE and 4.2 % 
for both the LSDA and ISOcc(c = 0.5). 

We define the averaged relative error in the ionization 
potential 


i5ip 


\ 


M 





(3) 


Here, the index j runs over all systems in the test set 
up to the total number M = 26, and e stands either for 
the shifted (ce-ho) or unshifted (£ho) highest occupied KS 
eigenvalue. 

Note that in Eq. Q the unsigned deviation from ex¬ 
perimental IPs is employed to avoid a misleading result 
of zero average relative error, which emerges when there 
occurs an overestimation for some systems and an un¬ 
derestimation for others. However, in order to be able to 
distinguish between systematic over- or underestimation, 
an additional measure is defined accordingly: 


1 ^ 

-^-p) • (4) 

While (5ip provides the mean deviation from experimen¬ 
tal values in %, the quantity S indicates the average 
trend of the prediction, being naturally confined to the 
interval [—1,1]. Namely, for a systematic overestimation 
we obtain S' = 1, and for a systematic underestimation 
S = —1. Both quantities, (5ip and S, were obtained 
for various DEAs and their ensemble-generalized coun¬ 
terparts. 

Fig. [3| shows the corresponding results for the LSDA, 
the semi-local PBE and BLYP, the global hybrid func¬ 
tionals B3LYP and PBEh(0.25), the EXX, and the 
ISOcc(0.5) local hybrid functional. Note that for EXX 
the results for the regular and ensemble-generalized func¬ 
tional coincide, because the Hartree-PEXX functional ex¬ 
hibits a zero potential shift (s^j- |l4Cll [ Fig. [4| provides 
the corresponding results for the PBEh(a) global hybrid 
functional as a function of the param eter a, i.e., on vari¬ 
ous amounts of non-local EXX |l4l| . The value a = 0 in 


this figure reproduces the PBE result. Eig. [5| depicts (5ip 
and S obtained with the local hybrid ISOcc. The latter 
functional was developed using the so-called local mixing 
function, rather than a fixed mixing ratio of non-local 
and semi-local components. It contains a free parameter 
c, which implicitly determines the intrinsic amount of 
EXX included in the local hybrid. Higher c values corre¬ 
spond to a higher fraction of EXX being included. Con¬ 
sequently, for ISOcc both the quantities given by Eqs. ® 
and o are functions of this parameter c, in analogy to 
the global hybrid PBEh(a). 



FIG. 3: Average relative error liip in % (upper axis) and 
signum function S (lower axis) for the LSDA, PBE, BLYP, 
B3LYP, PBE(a = 0.25), ISOcc(c = 0.5) and pure EXX 
(black) as well as their ensemble-generalized versions (red). 



0 - 0 0.2 0.4 0.6 0.8 1 

—I— I-1-1-1- 1 -1-1-1- 1 -1-1-1- 1 -1-1-1- 1 -1-1-1- 1 — a 



FIG. 4: Average relative error dip (a) in % (npper axis) and 
signum function S{a) (lower axis) for PBEh(a) (black) and 
e-PBEh(a) (red) in dependency on the parameter a. 
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FIG. 5: Average relative error (5ip(c) in % (upper axis) and 
signum function 5(c) (lower axis) for ISOcc(c) (black) and 
e-ISOcc(c) (red) in dependency on the parameter c. 

In principle, both parameters a and c, in PBEh(a) and 
ISOcc(c), respectively, are free. However, as mentioned 
earlier, it is known that in terms of total energy related 
quantities PBEh(a) performs best for a = 0.25, while 
we recently showed that for the ISOcc(c) functional the 
optimal parameter value is c = 0.5. Therefore, both func¬ 
tionals, using their optimal respective parameters play a 
special role in the following discussion, and Fig. [3] shows 
their performance in comparison to the other DFAs. 

Fig. [3] clearly indicates that using the unshifted eigen¬ 
values £ho, the three (semi-)local functionals LSDA, 
PBE, and BLYP strongly and systematically underes¬ 
timate i5ip by Ri 41 — 43%. Regarding hybrids, for the 
global hybrid B3LYP we obtain an underestimation of 
31%, for PBEh(a = 0.25) 28% and for the local hybrid 
ISOcc(c = 0.5) 26%. The improvement of hybrids over 
(semi-)local functionals is explained by the fact that the 
non-local terms in hybrids lead to a partial cancellation 
of the self-interaction error and an improved behavior of 
the xc potential in the asymptotic limit. 

If the parameters a and c are varied, Figs. [4] and [5] il¬ 
lustrate that when using £ho the global hybrid PBEh(a) 
and the local hybrid ISOcc(c) show a transition in their 
parameter-dependent S'-function from negative to posi¬ 
tive values. This feature clearly indicates that, for the 
systems studied here, it is possible to fit the correspond¬ 
ing functional parameter for a given system so that £ho 
exactly gives the experimental IP. If a and c are optimized 
to reduce the error dip, we obtain an underestimation of 
5 % for a = 0.75 in PBEh and of 6 % for c = 4.5 in ISOcc. 
Therefore, by changing the parameters a and c, we are 
able to strongly reduce the average error in the IP of our 
test set. However, this comes at a price in total energy- 
related quantities, as has been shown in Refs. (tI, [s^, 
and is the subject of the so-called ‘parameter dilemma’ 
presented in Sec. U 


When using the ensemble-corrected highest occupied 
eigenvalues Ee-hoj we obtain a completely different pic¬ 
ture. First, the systematic underestimation now changes 
to an overestimation. All of the aforementioned func¬ 
tionals now show a very similar average error of dip ~ 
14 — 17% , which is significantly smaller than the re¬ 
sults from non-generalized functionals. Second, for 
ensemble-generalized hybrid functionals e-PBEh(a) and 
e-ISOcc(c), there is no transition from an underestima¬ 
tion to overestimation regime, but rather a systematic 
overestimation of the IP, independent of the parameter 
value. In other words, the amount of non-locality in¬ 
cluded in the hybrid functional plays a minor role in the 
description of IPs via shifted KS eigenvalues, in contrast 
to their unshifted counterparts. This confirms that the 
mechanism of cancellation between the change in the po¬ 
tential shift of Eq. (12) and the highest occupied eigen¬ 
value with a varying amount of non-locality is not par¬ 
ticular to the O 2 molecule, but rather a systematic fea¬ 
ture of ensemble-generalized functionals. Furthermore, 
in the ensemble-generalized version of DFAs the ‘pa¬ 
rameter dilemma’ does not emerge: since the ensemble- 
generalized eigenvalues describe IPs with an accuracy al¬ 
most independent of the amount of EXX included, one 
cannot deduce a preferred value of the parameter by 
minimizing dip. Therefore, in principle one could use 
the functional with the parameter optimized to describe 
binding processes and structural quantities, and rely on 
the description of IPs via the shifted eigenvalues Ee-ho- 
In this case, our results for dip using thermochemically 
optimized functionals with ensemble-generalization (such 
as e-B3LYP, e-PBEh(a = 0.25) and e-ISOcc(c = 0.5)) 
indicate a clear improvement over their non-generalized 
counterparts. 

Our results further indicate that even functionals 
whose xc terms were constructed on different grounds 
and from different perspectives, such as for example the 
PBE and BLYP functional, yield similar values of roughly 
dip ~ 15 % after applying the ensemble generalization. 
As even the inclusion of non-local components does not 
lead to significant change, one might wonder if this ’’nat¬ 
ural border” of 15 % is inherent to the ensemble shift 
mechanism regardless of the specific form of the respec¬ 
tive DFA put to task. This question has been checked by 
varying the parameters p and k used in th e construction 
of the PBE exchange functional [13, Il09l| . We find that 
for different choices of n and k one obtains different val¬ 
ues for the average relative error dip. For instance, using 
PBE exchange with a value of /i = 1.0 together with the 
original k = 0.8401 results in an error of dip = 20% 
when using £e-ho, while a combination of the original 

= 0.21951 and n = 5.0 leads to dip = 8%. From 
this we conclude that the ensemble-generalization as such 
does not lead to a fixed systematic error in the descrip¬ 
tion of experimental IPs via KS eigenvalues. However, 
the results of this subsection suggest that after the en¬ 
semble generalization the functionals examined here have 
a common missing part, which causes the described dis- 








































































crepancy in (5ip. 

Before concluding this sub-section, we note that while 
we have focused our work on the IP of neutral atoms 
and diatomic molecules, IPs of ions may in principle be 
assessed in the same manner. In particular, the elec¬ 
tron affinity (EA) of the neutral can be explored as the 
IP of the singly charged anion (barring geometrical re¬ 
laxation). Unfortunately, for the atoms and very small 
molecules studied here, it is well-known [iiMa that 
with common semi-local approximations negative ions of 
small systems may erroneously be predicted to be un¬ 
stable. However, when pe rform ing calculations with fi¬ 
nite basis sets, as in, e.g., [ll3j| . unbound states can be 
artificially stabilized [114| . because the basis set effec¬ 
tively confines the unbound electron to the vicinity of 
the neutral system. Because the ensemble generalization 
discussed here does not change the total energies of sys¬ 
tems with integer electrons (including neutrals, cations, 
and anions), anions that are not bound with the un¬ 
derlying xc functional will remain unbound even if its 
ensemble-generalized version is employed. Furthermore, 
although ensemble-generalization will generally shift the 
energy levels, including the unoccupied ones, the ques¬ 
tion of whether the lowest unoccupied KS orbital has a 
bound or unbound character will not be affected [^ . 
because orbitals are unchanged by a uniform shift of the 
potential. 


C. Ensemble-generalization and the Aufbau 
principle 

In general, at zero temperature the energy levels in 
the KS system have to be occupied according to the Auf¬ 
bau principle, i.e., the levels are occupied without ’holes’, 
starting with the lowest ones up. In the following, we 
term such an occupation proper. An example for a proper 
occupation is given in Fig. [1] for the O 2 molecule. All 
calculations performed for this work, except for a few 
discussed below, yield proper occupation. 

In spin-polarized calculations a special situation can 
occur, when each of the spin channels is occupied prop¬ 
erly itself, while the system as a whole possesses a ’hole’ 
in its occupation. For example, this happens when the 
In level of the 4--channel appears lower than the ho level 
of the '('-channel. An occupation of this kind is termed 
proper in a broad sense. It is emphasized here that a 
broad-sense-proper density obeys all the required restric¬ 
tions related to a rigor ous defini tion and differentiability 
of energy functionals [HI [13 , therefore it can serve 
as a legitimate solution of a many-electron system. In 
the past, broad-sense-proper occupations have been ob¬ 
served in certain transition-metal and lanthanide atoms 
and ions in LSDA and PBE calcu lations [3 - as well as 
in the Li atom with the EXX [n3 and with the exact KS 
potential, which has been obtained fro m ac curate wave- 
function-methods based spin densities |ll8j |. The latter 
result strengthens our understanding that a broad-sense- 


proper occupation is not necessarily an artifact of some 
DFAs, but rather is an expected result, because it may 
appear even with the exact functional. 


-0- 

0 fs ry t 


- - 0 - _ 
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' ^ £e-lu 


up channel down channel 


FIG. 6: Diagram of the highest occupied and lowest unoccu¬ 
pied KS eigenvalues of the Li atom, for both spin channels, 
before and after application of the potential shifts of Eq. ©, 
obtained within the PBE functional, along with the negative 
of the experimental IP. All values are in Ha. The highest oc¬ 
cupied eigenvalue at the ^channel is lower than -0.65 Ha, and 
is therefore not shown for clarity. 

In the current work we find that the ensemble gen¬ 
eralization, by means of the potential shifts Vq, yields 
broad-sense-proper results for systems which appeared 
strictly proper before. Figure [6] illustrates the situation 
for the Li atom calculated with PBE and e-PBE. Due to 
the fact that Ug = —0.087 Ha, while Vq = —0.603 Ha, 
the ),-e-lu level appears below the f-e-ho level, causing 
a broad-sense-proper occupation. The significant differ¬ 
ences in the values of the two potential shifts are associ¬ 
ated with the different nature of the tr-ho orbitals: the 
f-ho orbital is a relatively delocalized, high-lying 2s or¬ 
bital, whereas the 4--ho orbital is a localized, low-lying Is 
orbital. A similar situation is observed in the Na atom, 
for which we obtained a broad-sense-proper result, too. 
To summarize, we view the appearance of broad-sense- 
proper occupations in the ensemble treatment of the al¬ 
kaline atoms Li and Na as another feature of the exact 
DFT result, which has been recovered by the ensemble 
generalization. 


D. The derivative discontinuity in 
ensemble-generalized functionals 

While in the current work we are concerned primarily 
with the IP of atoms and molecules, it is worth discussing 
a related quantity - the (fundamental) gap. Eg. By def¬ 
inition, Eg = I — A, where A is the electron affinity, i.e., 
the energy gained by adding one electron to the system. 
As opposed to the IP, the gap of the interacting system 
does not equal the gap of the KS system, E^ = £iu —£ho, 
even for the exact xc functional. Instead, Eg = E^ -b A, 
where A is the derivative discontinuity (DD) - a ’’jump” 
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experienced by the KS potential when it is varied with 
resp ect to N, and N crosses an integer value 

There exist several ways to find the DD. First, for finite 
systems it can be obtained using total energy differences: 


Ae = E{No + 1) - 2E{No) + E{No - 1) - Ef. (5) 
Seco nd, the DD can be obtained as suggested in [M 


AoEP — (^lu |'^xc,lu|^lu) (^lul'^xcl^lu) 


( 6 ) 


where Uxcy(r) := be., the orbital- 

specific xc potential of the i-th orbital and v^c ■= 
SExc/Sn^ i.e., the local xc potential, whi ch in gen eral has 
to be obtained via the OEP procedure |l0Cll - ll02| (hence 
the index OEP). The derivation of Aqep assumes the 
’’alignment equality” (iy9ho|Mxc.ho|</?ho) = (<Pho|fxc|</5ho), 
which determines the free constant in v-x_c(j) as part of 
the OEP procedure. Additional approache s to intro duce 
the derivative d iscon tinuity include Refs. [l28l - ll32 | (see 
also (dll. for an overview). 

Finally, an approximation for the DD can be obtained 
from an ensemble treatment for a given underlying Hxc 
functional, as proposed in Ref. [s^: 

Aens = Anxcln -b |</5iun — 2i?Hxc[n] -b i?Hxc[n — |:phon+ 
+ [ (fir WHxc[n](r) (|:/Jho(r)P - |:b’iu(r)P) (7) 


The first way requires three independent self-consistent 
calculations of the total energy (hence the index E): of the 
neutral system, the cation and the anion. In contrast, the 
second and third ways yield the DD from KS quantities 
of the neutral system only, which is an advantage when 
considering infinite systems. 

Relying on our experience with ensemble-generalized 
calculations for atoms and small molecules (Ref. [s^ and 
this work), we expect Aens obtained with an approximate 
xc functional to be larger than Ae. As has been shown 
in Fig. 2 of Ref. [s^ (lower panel), ffe-ho is obtained as 
being somewhat too low immediately to the left of an in¬ 
teger N and somewhat too high immediately to the right 
of it. As a result, Aens overestimates the true discontinu¬ 
ity. This overestimate is related to the residual concavity 
of the E{N) curve after the ensemble generalization. In 
the current study, we showed that the overestimate in 
Ee-ho to the left of the integer point, which corresponds 
to the negative of the IP, is systematic, i.e., it happens 
in various systems and with different functionals. Con¬ 
sequently, we expect a systematic overestimate for Agns 
and the resulting Eg. 

The discrepancy between Aqep and Agns has a differ¬ 
ent origin. While Aqep originates because the KS po¬ 
tentials are differently ’’aligned” (see above) to the left 
and to the right of an integer point, Agns comes from 
two sources (see Ref. [S^ for detailed explanations): the 
first is the same as for Aqep; the second is the fact that 


the ensemble-generalized KS potential does not approach 
zero at r —>■ oo, but rather a constant vq (see Eq. 
which is different to the left and to the right of an in¬ 
teger N. Aqep does not consider the second source 
descri bed abov e, assuming (correctly in the context of 
Refs. |l26l . Il27j l that the potentials asymptotically tend 
to zero. In fact, it can be analytically shown that Aqep 
is an ingredient in Ag„s, which was denoted by Ai in 
Ref. [13. 


V. CONCLUSIONS AND SUMMARY 


In the current work, we employed the ensemble- 
generalization procedure [13 for a test set of 26 diatomic 
molecules and single atoms, for a variety of xc func¬ 
tionals. These include the local spin-density approx¬ 
imation (LSDA), the semi-local PBE and BLYP, the 
global hybrids B3LYP and PBEh(a), and the local hy¬ 
brid ISOcc(c). We focused on the prediction of the IP 
via the highest occupied KS eigenvalue, Ee-ho- 

We found that implementing the ensemble approach 
improves, on average, the correspondence of Eho with the 
experimental IP for all xc functionals considered, chang¬ 
ing the general tendency in the IP prediction from a gross 
underestimation to a smaller overestimation. 

For functionals that include a parameter, namely the 
hybrids PBEh(a) and ISOcc(c), we observed a rather 
weak dependence of Eg-ho on the respective functional 
parameter, while yielding a roughly constant overestima¬ 
tion to the IP, with respect to experiment. This eases the 
so-called ‘parameter dilemma’: there are no two optimal 
values of the functional’s parameter originating from fit¬ 
ting to total energy-related quantities as opposed to fit¬ 
ting potential-related quantities. Instead, the parameter 
can be determined relying on energetics only, because of 
its weak influence on the value of Eg_ho- Indeed, the aver¬ 
age relative error in the ionization potential, 5ip, equals 
approximately 15 % for all ensemble-generalized xc func¬ 
tionals, as can be seen from Fig. [3l Surprisingly, such 
features of the underlying xc functional, as being local 
(LSDA), semi-local (PBE, BLYP) or non-local (B3LYP, 
PBEh, ISOcc), relying on features of the homogeneous 
electron gas (LSDA, PBE, PBEh, ISOcc) or not (BLYP, 
B3LYP), are of little relevance with respect to the IP 
prediction, once the functional is used in the ensemble- 
generalized form. We therefore conclude that upon en¬ 
semble generalization (Eq. ([T])) all the functionals we 
tested share the same deficiency. It is most probably 
related to the remaining concavity of the E{N) curve, 
due to the implicit dependence of the KS orbitals on a. 
Therefore, future improvement in the IP prediction via 
Eg-ho may be achieved via formulating a correction that 
will remove the remaining concavity in E(N). 
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arXiv:1508.07795vl [physics.chem-ph] 31 Aug 2015 


Supplemental Material to: “Effect of ensemble generalization on the 
highest-occupied Kohn-Sham eigenvalue” 

Eli Kraisler,^’El Tobias Schmidt,^ ’0 Stephan Kiimmel,^ and Leeor Kronik^ 

^Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100,Israel 
^Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany 

(Dated: July 5, 2017) 


In the following we present in tabular form the numerical data used to generate Figs. (2)-(5) of the main text. 
For each system in the test set, we provide the negative of the unshifted highest occupied eigenvalue, —Sho, as 
well as its shifted counterpart, —Se-ho, for a variety of functionals discussed in the article. Table U contains results 
obtained with the functionals LSDA, BLYP, B3LYP and pure EXX, tables HIl and IIIII provide values obtained with 
the global hybrid functional PBEh(a) and the local hybrid functional ISOcc(c), as function of the parameter a or c, 
respectively. Additionally, if calculated, the ionization potential evaluated by the ASCF method is included in the 
tables. Finally, for each functional the averaged relative errors d/p, evaluated with —Sho, —£e-ho and ASCF, together 
with the corresponding sign functions S, are presented. Our calculations were based on bond lengths and vertical 
ionization potentials determined by experiment as given in Ref. [ij and http: //webbook. nist. gov; the actual values 
are included in the tables for comparison. 


TABLE I: Unshifted and shifted highest occupied eigenvalue —Sho and —Se-ho for the system set evaluated with LSDA, BLYP, 
B3LYP and pure EXX. Additionally, for LSDA the IP computed via ASCF is listed. All energy values are in Ha . 


System 

Rab (Bohr) 

I Pexp 

IPdft via ... 

LSDA 

Functional 

BLYP B3LYP 

EXX 

Ha 

1.4011 

0.5669 

^ho 

0.3772 

0.3819 

0.4313 

0.5945 




^e—ho 

0.6257 

0.6329 

0.6353 

0.5945 




ASCF 

0.5963 




LiH 

3.0139 

0.285 

^ho 

0.1613 

0.1588 

0.1922 

0.3011 




^e — ho 

0.3424 

0.3435 

0.3437 

0.3011 




ASCF 

0.3017 




Lia 

5.0518 

0.1879 

£ho 

0.1189 

0.1127 

0.1313 

0.1812 




^e—ho 

0.2147 

0.2115 

0.2123 

0.1812 




ASCF 

0.1954 




LiF 

2.9553 

0.4155 

£ho 

0.2333 

0.2254 

0.2809 

0.4760 




^e—ho 

0.5835 

0.5733 

0.5622 

0.4760 




ASCF 

0.4542 




BeH 

2.5368 

0.3015 

£ho 

0.1692 

0.1666 

0.1983 

0.3096 




^e — ho 

0.3152 

0.3175 

0.3193 

0.3096 




ASCF 

0.3057 




BH 

2.3290 

0.359 

^ho 

0.2031 

0.2011 

0.2370 

0.3461 




^e—ho 

0.3876 

0.3858 

0.3869 

0.3461 




ASCF 

0.3558 




BO 

2.2766 

0.489 

^ho 

0.3045 

0.2462 

0.2857 

0.5199 




^e—ho 

0.5064 

0.4395 

0.4429 

0.5199 




ASCF 

0.4765 




BE 

2.3861 

0.4087 

£ho 

0.2508 

0.3004 

0.3491 

0.4053 




^e—ho 

0.4435 

0.5038 

0.5129 

0.4053 




ASCF 

0.4059 




CH 

2.1163 

0.391 

Sho 

0.2185 

0.2091 

0.2562 

0.4162 




Se—ho 

0.4520 

0.4411 

0.4432 

0.4162 




ASCF 

0.4057 




CN 

2.2144 

0.4997 

Sho 

0.3508 

0.3404 

0.3836 

0.5316 




Se—ho 

0.5932 

0.5822 

0.5778 

0.5316 




ASCF 

0.5376 




CO 

2.1322 

0.515 

Sho 

0.3350 

0.3311 

0.3839 

0.5526 
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TABLE I — continued from previous page 


System 

(Bohr) 

TP 

J J exp 

IPdft via ... 

LSDA 

Functional 

BLYP B3LYP 

EXX 




^e — ho 

0.5769 

0.5705 

0.5780 

0.5526 




ASCF 

0.5176 




NH 

1.9600 

0.4958 

^ho 

0.2928 

0.2842 

0.3404 

0.5011 




^e—ho 

0.5695^ 

0.5774^ 

0.573lt 

0.5011 




ASCF 

0.5066 




Nz 

2.0743 

0.5733 

^ho 

0.3825 

0.3770 

0.4364 

0.6303 




^e—ho 

0.6251 

0.6168 

0.6294 

0.6303 




ASCF 

0.5740 




NO 

2.1746 

0.3405 

^ho 

0.1690 

0.1622 

0.2206 

0.4173 




^e—ho 

0.4206 

0.4116 

0.4217 

0.4173 




ASCF 

0.3685 




OH 

1.8324 

0.4784 

£ho 

0.2740 

0.2629 

0.3203 

0.4984 




^e—ho 

0.5879 

0.5740 

0.5715 

0.4984 




ASCF 

0.4957 




O 2 

2.2819 

0.4531 

£ho 

0.2547 

0.2502 

0.3173 

0.5571 




^e—ho 

0.5316 

0.5265 

0.5396 

0.5571 




ASCF 

0.4683 




FH 

1.7326 

0.5924 

^ho 

0.3608 

0.3544 

0.4203 

0.6453 




^e—ho 

0.7391 

0.7313 

0.7235 

0.6453 




ASCF 

0.6172 




F 2 

2.6695 

0.5769 

^ho 

0.3544 

0.3481 

0.4205 

0.6667 




^e—ho 

0.6478 

0.6406 

0.6552 

0.6667 




ASCF 

0.5744 




H 


0.4997 

£ho 

0.2690 

0.2722 

0.3191 

0.5000 




^e—ho 

0.4787 

0.4979 

0.4991 

0.5000 




ASCF 

0.4787 




Li 


0.1981 

£ho 

0.1163 

0.1114 

0.1311 

0.1962 




^e—ho 

0.2013 

0.2034 

0.2041 

0.1962 




ASCF 

0.2011 




Be 


0.3426 

^ho 

0.2057 

0.2009 

0.2290 

0.3089 




^e—ho 

0.3447 

0.3439 

0.3454 

0.3089 




ASCF 

0.3318 




B 


0.3049 

^ho 

0.1509 

0.1494 

0.1866 

0.3170 




^e—ho 

0.3380 

0.3382 

0.3390 

0.3170 




ASCF 

0.3175 




C 


0.4138 

£ho 

0.2249 

0.2182 

0.2662 

0.4378 




^e—ho 

0.4714 

0.4667 

0.4661 

0.4378 




ASCF 

0.4313 




N 


0.5341 

£ho 

0.3085 

0.2970 

0.3560 

0.5705 




^e—ho 

0.6115 

0.5996 

0.5957 

0.5705 




ASCF 

0.5512 




0 


0.5005 

^ho 

0.2737 

0.2803 

0.3366 

0.5193 




^e—ho 

0.5962 

0.6045 

0.5984 

0.5193 




ASCF 

0.5146 




F 


0.6403 

^ho 

0.3808 

0.3797 

0.4478 

0.6779 




^e—ho 

0.7679 

0.7677 

0.7598 

0.6779 




ASCF 

0.6598 







Sip (%) via —Sho 

41.45 

42.62 

31.50 

9.19 




Sip (%) via -e^-ho 

16.42 

15.48 

15.48 

9.19 




Sip (%) via ASCF 

4.17 







S via —£ho 

- 1.00 

- 1.00 

- 1.00 

0.62 




S via Sg —fiQ 

0.92 

0.92 

0.92 

0.62 


t 


In this case, the states marking the highest occupied eigenvalue before and after the shift are not the same, but they belong to different 
spin channels. 



TABLE II: Unshifted and shifted highest occupied eigenvalue —Sho and —Ee-ho for the system set evaluated with the global hybrid functional PBEh(a) in dependence on 
the functional parameter a. Additionally, for pure PBE (i.e., a = 0) the IP computed via ASCF is listed. All energy values are in Ha . 


Functional parameter a 


System 

R7S (Bohr) 

TP 

exp 

IPdft via ... 

0 

0.1 

0.2 

0.25 

0.3 

0.4 

0.5 

0.6 

0.7 

0.75 

0.8 

0.9 

1.0 

H2 

1.4011 

0.5669 

^ho 

0.3815 

0.4053 

0.4292 

0.4412 

0.4532 

0.4772 

0.5013 

0.5256 

0.5498 

0.5620 

0.5742 

0.5986 

0.6230 




^e — ho 

0.6268 

0.6269 

0.6271 

0.6272 

0.6273 

0.6275 

0.6278 

0.6281 

0.6284 

0.6285 

0.6287 

0.6290 

0.6294 




ASCF 

0.5965 













LiH 

3.0139 

0.285 


0.1603 

0.1759 

0.1917 

0.1998 

0.2079 

0.2242 

0.2407 

0.2575 

0.2744 

0.2830 

0.2916 

0.3088 

0.3263 




^e — ho 

0.3393 

0.3382 

0.3372 

0.3367 

0.3363 

0.3355 

0.3347 

0.3340 

0.3333 

0.3330 

0.3327 

0.3321 

0.3316 




ASCF 

0.2959 













Li2 

5.0518 

0.1879 

^ho 

0.1185 

0.1268 

0.1351 

0.1392 

0.1434 

0.1518 

0.1602 

0.1687 

0.1772 

0.1814 

0.1857 

0.1943 

0.2029 




^e — ho 

0.2132 

0.2124 

0.2117 

0.2114 

0.2110 

0.2104 

0.2097 

0.2091 

0.2085 

0.2082 

0.2079 

0.2073 

0.2068 




ASCF 

0.1931 













LiF 

2.9553 

0.4155 

^ho 

0.2252 

0.2515 

0.2783 

0.2920 

0.3058 

0.3338 

0.3622 

0.3912 

0.4206 

0.4355 

0.4504 

0.4806 

0.5112 




^e — ho 

0.5731 

0.5660 

0.5592 

0.5559 

0.5527 

0.5464 

0.5403 

0.5344 

0.5287 

0.5259 

0.5232 

0.5179 

0.5127 




ASCF 

0.4437 













BeH 

2.5368 

0.3015 

^ho 

0.1728 

0.1876 

0.2024 

0.2098 

0.2173 

0.2322 

0.2472 

0.2622 

0.2773 

0.2848 

0.2924 

0.3075 

0.3227 




^e — ho 

0.3197 

0.3198 

0.3199 

0.3199 

0.3200 

0.3201 

0.3202 

0.3203 

0.3204 

0.3205 

0.3205 

0.3206 

0.3207 




ASCF 

0.3099 













BH 

2.3290 

0.359 

£ho 

0.2033 

0.2201 

0.2370 

0.2455 

0.2540 

0.2710 

0.2881 

0.3052 

0.3224 

0.3310 

0.3396 

0.3569 

0.3742 




^e — ho 

0.3846 

0.3839 

0.3833 

0.3829 

0.3826 

0.3820 

0.3813 

0.3807 

0.3801 

0.3798 

0.3796 

0.3790 

0.3785 




ASCF 

0.3520 













BO 

2.2766 

0.489 


0.2492 

0.2678 

0.2864 

0.2957 

0.3051 

0.3238 

0.3424 

0.3611 

0.3798 

0.3892 

0.3986 

0.4173 

0.4361 




^e — ho 

0.4382 

0.4386 

0.4389 

0.4391 

0.4392 

0.4395 

0.4398 

0.4400 

0.4402 

0.4404 

0.4404 

0.4406 

0.4408 




ASCF 

0.3997 













BE 

2.3861 

0.4087 

£ho 

0.3063 

0.3296 

0.3529 

0.3646 

0.3762 

0.3995 

0.4227 

0.4460 

0.4692 

0.4808 

0.4924 

0.5156 

0.5387 




^e — ho 

0.5063 

0.5100 

0.5135 

0.5153 

0.5170 

0.5203 

0.5235 

0.5266 

0.5296 

0.5311 

0.5325 

0.5353 

0.5380 




ASCF 

0.4770 













CH 

2.1163 

0.391 

^ho 

0.2155 

0.2378 

0.2602 

0.2715 

0.2828 

0.3056 

0.3284 

0.3514 

0.3745 

0.3860 

0.3976 

0.4208 

0.4442 




^e — ho 

0.4465 

0.4462 

0.4460 

0.4459 

0.4459 

0.4457 

0.4456 

0.4455 

0.4455 

0.4454 

0.4454 

0.4454 

0.4454 




ASCF 

0.4005 













CN 

2.2144 

0.4997 

^ho 

0.3452 

0.3655 

0.3859 

0.3961 

0.4064 

0.4271 

0.4480 

0.4694 

0.4912 

0.5023 

0.5136 

0.5368 

0.5608 




^e — ho 

0.5868 

0.5832 

0.5798 

0.5781 

0.5765 

0.5734 

0.5705 

0.5680 

0.5659 

0.5650 

0.5642 

0.5629 

0.5621 




ASCF 

0.5305 













CO 

2.1322 

0.515 

^ho 

0.3322 

0.3574 

0.3827 

0.3953 

0.4080 

0.4333 

0.4587 

0.4841 

0.5095 

0.5222 

0.5349 

0.5604 

0.5859 




^e — ho 

0.5685 

0.5709 

0.5732 

0.5744 

0.5755 

0.5777 

0.5799 

0.5821 

0.5842 

0.5852 

0.5862 

0.5882 

0.5902 




ASCF 

0.5092 













NH 

1.9600 

0.4958 


0.2911 

0.3179 

0.3449 

0.3584 

0.3720 

0.3994 

0.4270 

0.4518 

0.4744 

0.4858 

0.4971 

0.5199 

0.5428 




^e — ho 

0.5742^ 

0.5712t 

0.5683+ 

0.5669+ 

0.5656+ 

0.5629+ 

0.5602+ 

0.5576 

0.5551 

0.5539 

0.5528 

0.5504 

0.5481 




ASCF 

0.5079 













N 2 

2.0743 

0.5733 

£ho 

0.3773 

0.4058 

0.4343 

0.4486 

0.4629 

0.4916 

0.5204 

0.5491 

0.5780 

0.5924 

0.6068 

0.6357 

0.6647 




^e — ho 

0.6172 

0.6220 

0.6269 

0.6293 

0.6318 

0.6367 

0.6417 

0.6467 

0.6517 

0.6542 

0.6568 

0.6618 

0.6669 




ASCF 

0.5655 













NO 

2.1746 

0.3405 


0.1635 

0.1914 

0.2195 

0.2337 

0.2479 

0.2764 

0.3050 

0.3338 

0.3627 

0.3771 

0.3916 

0.4207 

0.4498 




^e — ho 

0.4125 

0.4162 

0.4200 

0.4219 

0.4238 

0.4276 

0.4314 

0.4353 

0.4392 

0.4411 

0.4431 

0.4469 

0.4508 




ASCF 

0.3614 
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TABLE II — continued from previous page 


System 

(Bohr) 

TP 

J- exp 

IPdft via ... 

0 

0.1 

0.2 

0.25 

Functional parameter a 

0.3 0.4 0.5 0.6 

0.7 

0.75 

0.8 

0.9 

1.0 

OH 

1.8324 

0.4784 

^ho 

0.2636 

0.2908 

0.3182 

0.3320 

0.3458 

0.3736 

0.4016 

0.4297 

0.4579 

0.4721 

0.4863 

0.5148 

0.5435 




^e — ho 

0.5727 

0.5697 

0.5668 

0.5654 

0.5639 

0.5612 

0.5585 

0.5558 

0.5533 

0.5520 

0.5508 

0.5483 

0.5460 




ASCF 

0.4804 













O 2 

2.2819 

0.4531 

£ho 

0.2507 

0.2829 

0.3155 

0.3319 

0.3484 

0.3815 

0.4150 

0.4486 

0.4825 

0.4995 

0.5166 

0.5509 

0.5854 




^e — ho 

0.5261 

0.5315 

0.5370 

0.5398 

0.5427 

0.5485 

0.5545 

0.5606 

0.5668 

0.5699 

0.5731 

0.5796 

0.5862 




ASCF 

0.4635 













FH 

1.7326 

0.5924 

^ho 

0.3547 

0.3864 

0.4183 

0.4344 

0.4505 

0.4830 

0.5157 

0.5487 

0.5818 

0.5985 

0.6152 

0.6487 

0.6825 




^e — ho 

0.7312 

0.7260 

0.7209 

0.7183 

0.7158 

0.7110 

0.7062 

0.7016 

0.6971 

0.6948 

0.6926 

0.6883 

0.6841 




ASCF 

0.6093 













F 2 

2.6695 

0.5769 

^ho 

0.3475 

0.3824 

0.4174 

0.4351 

0.4528 

0.4883 

0.5240 

0.5599 

0.5960 

0.6142 

0.6323 

0.6688 

0.7054 




^e — ho 

0.6398 

0.6457 

0.6518 

0.6549 

0.6580 

0.6644 

0.6711 

0.6778 

0.6847 

0.6882 

0.6917 

0.6989 

0.7062 




ASCF 

0.5669 













H 


0.4997 


0.2791 

0.3016 

0.3242 

0.3355 

0.3469 

0.3698 

0.3928 

0.4159 

0.4392 

0.4508 

0.4625 

0.4860 

0.5096 




^e — ho 

0.5000 

0.5005 

0.5010 

0.5013 

0.5016 

0.5022 

0.5028 

0.5034 

0.5040 

0.5043 

0.5047 

0.5053 

0.5060 




ASCF 

0.5000 













Li 


0.1981 


0.1186 

0.1275 

0.1363 

0.1407 

0.1451 

0.1540 

0.1628 

0.1716 

0.1803 

0.1847 

0.1891 

0.1979 

0.2066 




^e — ho 

0.2055 

0.2052 

0.2049 

0.2048 

0.2047 

0.2044 

0.2041 

0.2039 

0.2036 

0.2035 

0.2034 

0.2032 

0.2029 




ASCF 

0.2053 













Be 


0.3426 

^ho 

0.2061 

0.2190 

0.2320 

0.2384 

0.2449 

0.2578 

0.2708 

0.2837 

0.2966 

0.3031 

0.3095 

0.3224 

0.3353 




^e — ho 

0.3436 

0.3431 

0.3427 

0.3425 

0.3423 

0.3419 

0.3414 

0.3410 

0.3406 

0.3404 

0.3402 

0.3399 

0.3395 




ASCF 

0.3307 













B 


0.3049 

^ho 

0.1534 

0.1708 

0.1884 

0.1973 

0.2062 

0.2241 

0.2422 

0.2603 

0.2786 

0.2878 

0.2970 

0.3156 

0.3342 




^e — ho 

0.3392 

0.3387 

0.3383 

0.3381 

0.3378 

0.3374 

0.3371 

0.3368 

0.3365 

0.3363 

0.3362 

0.3360 

0.3358 




ASCF 

0.3186 













C 


0.4138 


0.2241 

0.2467 

0.2695 

0.2810 

0.2925 

0.3157 

0.3390 

0.3626 

0.3862 

0.3981 

0.4101 

0.4340 

0.4581 




^e — ho 

0.4697 

0.4683 

0.4670 

0.4663 

0.4657 

0.4645 

0.4633 

0.4622 

0.4612 

0.4607 

0.4602 

0.4592 

0.4583 




ASCF 

0.4308 













N 


0.5341 

£ho 

0.3052 

0.3330 

0.3611 

0.3752 

0.3894 

0.4179 

0.4467 

0.4756 

0.5048 

0.5195 

0.5342 

0.5637 

0.5934 




^e — ho 

0.6058 

0.6007 

0.5999 

0.5995 

0.5992 

0.5985 

0.5979 

0.5974 

0.5970 

0.5968 

0.5966 

0.5963 

0.5961 




ASCF 

0.5414 













0 


0.5005 

^ho 

0.2794 

0.3067 

0.3343 

0.3481 

0.3621 

0.3901 

0.4183 

0.4467 

0.4752 

0.4896 

0.5040 

0.5329 

0.5619 




^e — ho 

0.6014 

0.5974 

0.5936 

0.5917 

0.5899 

0.5862 

0.5827 

0.5792 

0.5759 

0.5743 

0.5726 

0.5695 

0.5664 




ASCF 

0.5167 













F 


0.6403 

^ho 

0.3788 

0.4116 

0.4448 

0.4614 

0.4782 

0.5118 

0.5458 

0.5799 

0.6142 

0.6315 

0.6489 

0.6837 

0.7186 




^e — ho 

0.7663 

0.7611 

0.7561 

0.7537 

0.7512 

0.7465 

0.7419 

0.7374 

0.7329 

0.7308 

0.7287 

0.7245 

0.7204 




ASCF 

0.6488 
















Sip (%) via -Sho 

41.72 

36.44 

31.13 

28.48 

25.82 

20.52 

15.30 

10.35 

6.35 

5.45 

5.87 

9.51 

14.47 




Sip (%) via -Se-ho 

15.48 

15.20 

14.99 

14.90 

14.83 

14.71 

14.63 

14.60 

14.61 

14.63 

14.66 

14.75 

14.88 




Sip (%) via ASCF 

3.41 
















S via -£ho 

- 1.00 

- 1.00 

- 1.00 

- 1.00 

- 1.00 

- 1.00 

- 1.00 

- 1.00 

- 0.62 

- 0.38 

0.31 

0.69 

0.92 




S via Se—ixo 

1.00 

1.00 

1.00 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 

0.92 


t 




In this case, the states marking the highest occupied eigenvalue before and after the shift are not the same, but they belong to different spin-channels. 



TABLE III: Unshifted and shifted highest occupied eigenvalue —Eho and —e^-ho for the system set evaluated with the local hybrid functional ISOcc(c) in dependence on 
the functional parameter c. Additionally, the IP computed via ASCF is listed. All energy values are in Ha . 


Functional parameter c 


System 

(Bohr) 

TP 

J J exp 

IPdft via ... 

0 

0.1 

0.3 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

H2 

1.4011 

0.5669 

^ho 

0.3772 

0.3954 

0.4249 

0.4475 

0.4860 

0.5108 

0.5282 

0.5413 

0.5515 

0.5597 

0.5665 

0.5722 

0.5771 




^e—ho 

0.6036 

0.6122 

0.6246 

0.6332 

0.6465 

0.6543 

0.6594 

0.6631 

0.6659 

0.6680 

0.6697 

0.6711 

0.6723 




ASCF 

0.5672 

0.5782 


0.6024 

0.6169 

0.6252 

0.6306 

0.6345 

0.6374 

0.6396 

0.6414 

0.6429 

0.6442 

LiH 

3.0139 

0.2850 

£ho 

0.1614 

0.1673 

0.1814 

0.1939 

0.2182 

0.2355 

0.2485 

0.2585 

0.2666 

0.2732 

0.2787 

0.2834 

0.2874 




^e—ho 

0.3377 

0.3325 

0.3359 

0.3394 

0.3462 

0.3508 

0.3541 

0.3566 

0.3585 

0.3601 

0.3614 

0.3624 

0.3633 




ASCF 

0.2814 

0.2829 


0.2940 

0.3032 

0.3092 

0.3134 

0.3166 

0.3191 

0.3212 

0.3229 

0.3244 

0.3257 

Li2 

5.0518 

0.1879 

^ho 

0.1189 

0.1202 

0.1246 

0.1286 

0.1368 

0.1430 

0.1480 

0.1522 

0.1557 

0.1588 

0.1614 

0.1638 

0.1659 




^e—ho 

0.2096 

0.2073 

0.2076 

0.2088 

0.2116 

0.2137 

0.2154 

0.2167 

0.2178 

0.2187 

0.2195 

0.2201 

0.2207 




ASCF 

0.1875 

0.1855 


0.1884 

0.1918 

0.1943 

0.1961 

0.1976 

0.1988 

0.1998 

0.2007 

0.2014 

0.2021 

LiF 

2.9553 

0.4155 

^ho 

0.2333 

0.2449 

0.2708 

0.2935 

0.3363 

0.3655 

0.3867 

0.4029 

0.4156 

0.4259 

0.4344 

0.4416 

0.4477 




^e—ho 

0.5808 

0.5636 

0.5511 

0.5463 

0.5423 

0.5412 

0.5407 

0.5405 

0.5402 

0.5400 

0.5398 

0.5396 

0.5395 




ASCF 

0.4518 

0.4355 


0.4200 

0.4173 

0.4170 

0.4172 

0.4176 

0.4180 

0.4184 

0.4189 

0.4193 

0.4198 

BeH 

2.5368 

0.3015 

£ho 

0.2441 

0.2471 

0.2543 

0.2603 

0.2711 

0.2782 

0.2833 

0.2871 

0.2901 

0.2925 

0.2945 

0.2963 

0.2977 




^e—ho 

0.3263 

0.3185 

0.3142 

0.3132 

0.3133 

0.3140 

0.3148 

0.3154 

0.3159 

0.3164 

0.3168 

0.3171 

0.3174 




ASCF 

0.3116 

0.3040 


0.3013 

0.3027 

0.3039 

0.3049 

0.3056 

0.3062 

0.3067 

0.3071 

0.3075 

0.3077 

BH 

2.3290 

0.359 

^ho 

0.2031 

0.2121 

0.2281 

0.2412 

0.2654 

0.2822 

0.2946 

0.3043 

0.3121 

0.3185 

0.3239 

0.3286 

0.3325 




^e—ho 

0.3818 

0.3802 

0.3823 

0.3851 

0.3910 

0.3949 

0.3978 

0.3999 

0.4015 

0.4028 

0.4039 

0.4047 

0.4055 




ASCF 

0.3477 

0.3478 


0.3536 

0.3594 

0.3632 

0.3659 

0.3679 

0.3693 

0.3705 

0.3715 

0.3723 

0.3729 

BO 

2.2766 

0.489 

£ho 

0.3499 

0.3592 

0.3786 

0.3953 

0.4265 

0.4482 

0.4599 

0.4687 

0.4756 

0.4812 

0.4858 

0.4897 

0.4930 




^e—ho 

0.5250^ 

0.5152+ 

0.5110+ 

0.5114+ 

0.5155+ 

0.5193+ 

0.5222+ 

0.5245+ 

0.5263+ 

0.5277+ 

0.5289+ 

0.5298+ 

0.5306+ 




ASCF 

0.4827 

0.4771 


0.4799 

0.4858 

0.4902 

0.4933 

0.4957 

0.4976 

0.4991 

0.5004 

0.5014 

0.5023 

BF 

2.3861 

0.4087 

^ho 

0.2509 

0.2605 

0.2778 

0.2920 

0.3183 

0.3364 

0.3498 

0.3601 

0.3684 

0.3753 

0.3810 

0.3859 

0.3901 




^e—ho 

0.4345 

0.4336 

0.4368 

0.4407 

0.4483 

0.4534 

0.4569 

0.4595 

0.4615 

0.4631 

0.4644 

0.4654 

0.4663 




ASCF 

0.3953 

0.3960 


0.4036 

0.4111 

0.4162 

0.4197 

0.4224 

0.4244 

0.4261 

0.4274 

0.4285 

0.4295 

CH 

2.1163 

0.391 

£ho 

0.2395 

0.2459 

0.2620 

0.2762 

0.3031 

0.3218 

0.3357 

0.3465 

0.3551 

0.3623 

0.3682 

0.3734 

0.3778 




^e—ho 

0.4518 

0.4386 

0.4332 

0.4327 

0.4342 

0.4360 

0.4375 

0.4387 

0.4396 

0.4404 

0.4410 

0.4416 

0.4421 




ASCF 

0.4099 

0.3978 


0.3916 

0.3925 

0.3939 

0.3952 

0.3962 

0.3971 

0.3979 

0.3985 

0.3990 

0.3995 

CN 

2.2144 

0.4997 

^ho 

0.3524 

0.3601 

0.3774 

0.3923 

0.4203 

0.4398 

0.4543 

0.4656 

0.4746 

0.4822 

0.4884 

0.4938 

0.4985 




^e—ho 

0.5942 

0.5841 

0.5786 

0.5771 

0.5766 

0.5770 

0.5773 

0.5775 

0.5776 

0.5777 

0.5777 

0.5777 

0.5777 




ASCF 

0.5433 

0.5416 


0.5521 

0.5634 

0.5714 

0.5773 

0.5820 

0.5857 

0.5887 

0.5913 

0.5934 

0.5953 

CO 

2.1322 

0.515 

^ho 

0.3350 

0.3496 

0.3745 

0.3946 

0.4309 

0.4556 

0.4737 

0.4876 

0.4988 

0.5079 

0.5155 

0.5219 

0.5275 




^e—ho 

0.5717 

0.5702 

0.5743 

0.5794 

0.5896 

0.5965 

0.6014 

0.6050 

0.6077 

0.6099 

0.6116 

0.6130 

0.6142 




ASCF 

0.5126 

0.5134 


0.5231 

0.5326 

0.5390 

0.5434 

0.5467 

0.5492 

0.5512 

0.5529 

0.5542 

0.5553 

NH 

1.9600 

0.4958 

£ho 

0.3157 

0.3320 

0.3573 

0.3770 

0.4093+ 

0.4305+ 

0.4461+ 

0.4581+ 

0.4678+ 

0.4757+ 

0.4824+ 

0.4880+ 

0.4929+ 




^e—ho 

0.5586 

0.5571 

0.5590 

0.5618 

0.5632 

0.5633 

0.5636 

0.5638 

0.5640 

0.5641 

0.5643 

0.5644 

0.5644 




ASCF 

0.4846 

0.4872 


0.4968 

0.5047 

0.5099 

0.5134 

0.5161 

0.5181 

0.5196 

0.5209 

0.5219 

0.5228 

N 2 

2.0743 

0.5733 

£ho 

0.3825 

0.3972 

0.4239 

0.4456 

0.4849 

0.5115 

0.5311 

0.5463 

0.5584 

0.5684 

0.5767 

0.5839 

0.5901 




^e—ho 

0.6252 

0.6229 

0.6271 

0.6327 

0.6443 

0.6526 

0.6586 

0.6631 

0.6667 

0.6696 

0.6719 

0.6739 

0.6755 




ASCF 

0.5745 

0.5735 


0.5846 

0.5964 

0.6046 

0.6105 

0.6149 

0.6184 

0.6211 

0.6234 

0.6253 

0.6270 

NO 

2.1746 

0.3405 

^ho 

0.1758 

0.1865 

0.2116 

0.2331 

0.2728 

0.2997 

0.3194 

0.3345 

0.3466 

0.3564 

0.3646 

0.3715 

0.3775 




^e—ho 

0.4184 

0.4092 

0.4093 

0.4132 

0.4228 

0.4298 

0.4349 

0.4388 

0.4418 

0.4442 

0.4462 

0.4478 

0.4492 




ASCF 

0.3698 

0.3615 


0.3659 

0.3752 

0.3819 

0.3868 

0.3905 

0.3934 

0.3958 

0.3977 

0.3993 

0.4007 


Continued on next page 



TABLE III — continued from previous page 


System 

(Bohr) 

TP 

J J exp 

IPdft via ... 

0 

0.1 

0.3 

0.5 

1.0 

Functional parameter c 
1.5 2.0 2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

OH 

1.8324 

0.4784 

^ho 

0.2747 

0.2911 

0.3208 

0.3448 

0.3875 

0.4158 

0.4362 

0.4516 

0.4638 

0.4736 

0.4818 

0.4887 

0.4946 




^e—ho 

0.5826 

0.5711 

0.5667 

0.5669 

0.5702 

0.5731 

0.5752 

0.5768 

0.5781 

0.5790 

0.5798 

0.5804 

0.5809 




ASCF 

0.4815 

0.4749 


0.4760 

0.4814 

0.4853 

0.4880 

0.4900 

0.4915 

0.4926 

0.4935 

0.4942 

0.4948 

O 2 

2.2819 

0.4531 

£ho 

0.2641 

0.2780 

0.3078 

0.3329 

0.3787 

0.4098 

0.4324 

0.4498 

0.4635 

0.4748 

0.4841 

0.4921 

0.4989 




^e—ho 

0.5300 

0.5219 

0.5232 

0.5281 

0.5398 

0.5484 

0.5547 

0.5595 

0.5633 

0.5663 

0.5688 

0.5709 

0.5727 




ASCF 

0.4697 

0.4628 


0.4699 

0.4816 

0.4902 

0.4965 

0.5014 

0.5052 

0.5083 

0.5109 

0.5130 

0.5149 

FH 

1.7326 

0.5924 

^ho 

0.3608 

0.3792 

0.4129 

0.4403 

0.4895 

0.5221 

0.5455 

0.5633 

0.5772 

0.5885 

0.5978 

0.6057 

0.6124 




^e—ho 

0.7363 

0.7222 

0.7135 

0.7108 

0.7095 

0.7096 

0.7098 

0.7100 

0.7101 

0.7101 

0.7101 

0.7101 

0.7101 




ASCF 

0.6179 

0.6250 


0.6576 

0.6845 

0.7024 

0.7152 

0.7248 

0.7322 

0.7382 

0.7430 

0.7471 

0.7505 

F 2 

2.6695 

0.5769 

^ho 

0.3544 

0.3741 

0.4118 

0.4425 

0.4972 

0.5333 

0.5591 

0.5785 

0.5937 

0.6060 

0.6162 

0.6247 

0.6320 




^e—ho 

0.6473 

0.6411 

0.6447 

0.6513 

0.6660 

0.6763 

0.6839 

0.6895 

0.6939 

0.6975 

0.7004 

0.7028 

0.7048 




ASCF 

0.5748 

0.5848 


0.6293 

0.6660 

0.6901 

0.7072 

0.7200 

0.7300 

0.7381 

0.7447 

0.7502 

0.7549 

H 


0.4997 

^ho 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 




^e—ho 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 




ASCF 

0.5000 

0.5000 


0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

0.5000 

Li 


0.1981 

£ho 

0.1788 

0.1780 

0.1786 

0.1797 

0.1821 

0.1838 

0.1852 

0.1863 

0.1873 

0.1881 

0.1888 

0.1894 

0.1900 




^e — ho 

0.2058 

0.1982 

0.1962 

0.1963 

0.1972 

0.1978 

0.1983 

0.1987 

0.1990 

0.1993 

0.1995 

0.1996 

0.1997 




ASCF 

0.2056 

0.1979 


0.1962 

0.1971 

0.1978 

0.1983 

0.1986 

0.1989 

0.1991 

0.1993 

0.1994 

0.1996 

Be 


0.3426 

^ho 

0.2056 

0.2118 

0.2225 

0.2312 

0.2473 

0.2586 

0.2671 

0.2739 

0.2794 

0.2841 

0.2881 

0.2915 

0.2945 




^e—ho 

0.3360 

0.3351 

0.3381 

0.3411 

0.3466 

0.3504 

0.3531 

0.3551 

0.3567 

0.3580 

0.3591 

0.3600 

0.3608 




ASCF 


0.3227 


0.3297 

0.3354 

0.3391 

0.3417 

0.3437 

0.3452 

0.3465 

0.3475 

0.3484 

0.3491 

B 


0.3049 

^ho 

0.1894 

0.1925 

0.2015 

0.2100 

0.2270 

0.2395 

0.2490 

0.2566 

0.2628 

0.2679 

0.2723 

0.2761 

0.2793 




^e—ho 

0.3377 

0.3277 

0.3236 

0.3234 

0.3252 

0.3270 

0.3285 

0.3296 

0.3305 

0.3312 

0.3318 

0.3323 

0.3327 




ASCF 

0.3213 

0.3116 


0.3067 

0.3078 

0.3093 

0.3105 

0.3115 

0.3122 

0.3128 

0.3133 

0.3137 

0.3141 

C 


0.4138 

£ho 

0.2740 

0.2802 

0.2943 

0.3067 

0.3302 

0.3468 

0.3592 

0.3688 

0.3766 

0.3830 

0.3884 

0.3930 

0.3970 




^e—ho 

0.4737 

0.4617 

0.4555 

0.4540 

0.4536 

0.4541 

0.4546 

0.4549 

0.4552 

0.4555 

0.4557 

0.4558 

0.4560 




ASCF 

0.4384 

0.4269 


0.4184 

0.4175 

0.4176 

0.4179 

0.4181 

0.4183 

0.4184 

0.4185 

0.4186 

0.4186 

N 


0.5341 

£ho 

0.3685 

0.3780 

0.3974 

0.4138 

0.4438 

0.4645 

0.4796 

0.4912 

0.5005 

0.5082 

0.5145 

0.5199 

0.5246 




^e—ho 

0.6158 

0.6019 

0.5944 

0.5923 

0.5913 

0.5922 

0.5928 

0.5933 

0.5937 

0.5940 

0.5942 

0.5944 

0.5946 




ASCF 

0.5590 

0.5448 


0.5356 

0.5349 

0.5353 

0.5358 

0.5362 

0.5365 

0.5367 

0.5369 

0.5371 

0.5372 

0 


0.5005 

^ho 

0.2753 

0.2941 

0.3272 

0.3536 

0.4001 

0.4308 

0.4526 

0.4690 

0.4819 

0.4924 

0.5010 

0.5083 

0.5145 




^e—ho 

0.5842 

0.5771 

0.5769 

0.5799 

0.5872 

0.5924 

0.5960 

0.5987 

0.6007 

0.6022 

0.6035 

0.6046 

0.6054 




ASCF 

0.4947 

0.4905 


0.4974 

0.5062 

0.5121 

0.5162 

0.5192 

0.5214 

0.5231 

0.5245 

0.5257 

0.5266 

F 


0.6403 

^ho 

0.3810 

0.4025 

0.4414 

0.4724 

0.5269 

0.5622 

0.5872 

0.6060 

0.6206 

0.6324 

0.6420 

0.6502 

0.6570 




^e—ho 

0.7626 

0.7491 

0.7428 

0.7424 

0.7449 

0.7473 

0.7490 

0.7502 

0.7511 

0.7518 

0.7523 

0.7528 

0.7531 




ASCF 

0.6464 

0.6362 


0.6358 

0.6415 

0.6457 

0.6485 

0.6505 

0.6520 

0.6531 

0.6539 

0.6546 

0.6551 




Sip (%) via -£ho 

36.50 

34.16 

29.56 

25.74 

18.89 

14.40 

11.31 

9.14 

7.67 

6.73 

6.21 

6.03 

6.07 




Sip (%) via -£e-ho 

15.83 

14.04 

13.56 

13.81 

14.73 

15.50 

16.09 

16.54 

16.89 

17.18 

17.41 

17.61 

17.78 




Sip (%) via ASCF 

4.13 

3.19 


4.23 

6.05 

7.39 

8.37 

9.13 

9.72 

10.20 

10.60 

10.94 

11.22 




S via -£ho 

- 0.92 

- 0.92 

- 0.92 

- 0.92 

- 0.92 

- 0.92 

- 0.92 

- 0.85 

- 0.62 

- 0.62 

- 0.15 

0.00 

0.08 




S via Se — ho 

0.92 

0.92 

0.85 

0.85 

0.92 

0.92 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


^ In this case, the states marking the highest occupied eigenvalue before and after the shift are not the same, but they belong to different spin channels. 

^ Here, the states marking the highest occupied eigenvalue are of a different spin channel compared to the results obtained with the parameter c = 0, i.e., a transition occurs. 
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